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We study the linear response to an external electric field of a system of fermions in a lattice at 
1 zero temperature. This allows to measure numerically the Euclidean conductivity which turns out to 

be compatible with an analytical calculation for free fermions. The numerical method is generalizable 
on : to systems with dynamical interactions where no analytical approach is possible. 

ON 

g . 

i ■ 1 Introduction 

The study of transport properties in a system of charged fermions is an interesting subject in areas 
as different as Physics of Plasma, Quantum Chromodynamics or Metals. In particular, the measure 
of the electrical conductivity is a very difficult and yet interesting problem, specially in presence of 
nonperturbative effects. In such a case, numerical methods are called for. In this work, we want to 
show that lattice-regularized Euclidean field theories can be useful in this respect, at least in the limit of 
vanishing temperature. However, it should be emphasized that the so called sign-problem needs still to 
be overcome in many cases (see, though, Refs. [| | for some successful simulations at finite density). 

In this paper, we restrict ourselves to a model consisting of fermions that only interact with an external 
electromagnetic field. In spite of its simplicity, it shares many properties with more realistic models, and 
it can be considered as a necessary first step to check any numerical method which could be used in 
models with dynamical interactions. We consider the standard U{\) lattice-action, with Wilson fermions 
and finite chemical potential, but with the gauge variables held fixed. We shall study the residue of the 
pole of the electrical conductivity at zero-frequency, which is purely imaginary. Since a non-vanishing 
value for this residue unambiguously signals a conducting phase, this is a rather interesting quantity in 
our opinion. In order to obtain it, we measure the electrical-current induced in the system by an external 
electric field. This technique requires a numerical calculation even in the case of an external spatially- 
homogeneous time-dependent electromagnetic field. The delicate point, however, is that our electric 
field varies in Euclidean time. One can nevertheless assume that there is a linear relation between the 
Euclidean current and the Euclidean electric field, at least for small fields. This Euclidean conductivity 
presents a pole whose residue can be straightforwardly measured. To check that the obtained result 
is physical, we follow a very elegant procedure due to KohnQ]. He showed that the real-time residue 
can be measured by studying the sensitivity of the ground-state energy to an external Aharonov-Bohm 
electromagnetic field. We show how can this be done in the lattice formalism, and, in this particularly 



'Electronic address 
T Electronic address 
t Electronic address 



buj@gteorico .unizar . es 
laf Olattice . f is .ucm. es 
Victor .Mart inOromal . inf n. it 



1 



simple case of free fermions, we calculate it (unfortunately, the Kohn recipe seems really hard to use in 
a Monte Carlo study of a self-interacting problem) . Although at present we lack a rigorous proof of the 
equivalence of both calculations, its excellent numerical agreement gives a strong support to the linear 
response method. 



2 The Model 

Let us consider a model of Wilson Fermions 0, ^| in a lattice of spacing a s in the three spatial directions 
and at in the temporal one, coupled to an external electromagnetic field. We denote by A to the quotient 
a t /a s . The partition function can be written as (the * superscript stands for complex conjugation) @, |8[ §| 

^2%M xy (U)% , (1) 

.x,y 



Z[U] 



J Yl d* z d* z exp 



M xy (U) = e^U x ,o(j -r t )5^ x+6 -e-^Ul ( l0 + r t )6 y+ - 0tX (2) 

3 

+ A^[[/ X ,i(7i - r s )5y iX+i - U*^(ji + r s )5 y+i . x ] + [(2m + 6r s )A + 2r t ]S x . y , 

i=l 

where U x ,u = e lAx - v , A being the gauge field, and ^ x are the anticonmuting Grassmann fermionic 
fields. The indices x, y run on the points of the four-dimensional space-time lattice. We impose periodic 
boundary conditions for the gauge field, and periodic in space but antiperiodic in time [v = 0) for the 
Grassmann field. The site x + v is the neighbor of x in the v = 0, 1, 2, 3 direction. For finite temporal 
length, Lq, the system is at finite temperature T — (a t £o) • In this paper we will only consider the zero 
temperature (Lq — ► oo) limit. We follow the prescription of introducing the chemical potential through 
an imaginary gauge field A = (— iA/i, 0, 0, 0) ||, which is fairly convenient for analytical calculations. 

The Wilson parameter, r, can be taken different for the spatial and time directions. In the limit A — > 
with a s fixed the model describes a spatial lattice with a continuous time (as electrons in a metal) , while 
for a continuum field theory both spatial and time continuum limits should be taken. 

We shall use the following representation for the (Euclidean) gamma matrices 
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where <Ji are the Pauli matrices. 

To define the electric four-current in the lattice we recall that in the space-time continuum limit it is 
defined as 

j v (x) =^(x)j v ^(x), (4) 

that can be obtained as a logarithmic derivative of the partition function respect to the gauge-field. This 
calculation can be exactly mimicked on the lattice noticing that a change in the link variable should be 
of the form U XM — » e lctx -"U Xt „. In this way one obtains ||: 

.SlogZ 

{jx.u) =1-^ , (5) 

oa x ^ v 

where now 

(j*,o) - (*,e^t/ :c ,o(7o-r t )* x+ e + * :c+ ae-^f/*, (7o+rt)* :c ) (6) 
(j x ,i) = A(* x f/x,i(7i -r^x+i + Vx+iU^iiji +r s )* x > , i = 1,2,3. (7) 

The jo component is just the electric charge density that one encounters by differentiating with respect to 
Xfi the free energy density S. Moreover, from the gauge invariance of the determinant of the fermionic 
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matrix, M, it is straightforward to prove the lattice continuity equation, for any configuration of the 
electromagnetic-field : 

- Yl ~ ■ (8) 



Eqs. (^|) and (Q) can be written free of Grassmann variables as 

- (j*,o) = e A ^,oTr[( 7 o - n)M^J + e^U^Tr^ + r t )M^ +0J 
- (jx,i> = \U X:i Tr{('fi - r s )M-^J + XU* ti Tr[( 7i + r s )M~^ , 
where Tr stands for the trace over Dirac indices. The above expressions and the relation 

M{U*) = 7173 {M{U)f 7371 , 

allow to prove that 

(jx,u)u = {jx,v)u* ■ 

In an uniform electrical field, the charge density should remain constant under field inversion, while the 
electrical current should change sign. Therefore, from Eq. (|l|) one expects the former to be real and the 
latter to be imaginary (Euclidean space-time!). 

In absence of external fields (U = 1) the matrix M can be diagonalized in Fourier space, which allows 
to explicitly perform the functional integrals, and to compute the free energy and the propagator. For 
brevity, we only quote the result for the charge density in the case ft = r s = 1, fi > 0, that in the infinite 
volume limit reads (see Refs. jl0[ |llj for similar calculations), 



(9) 
(10) 

(11) 
(12) 



where 
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d J fc 



6(n- X- 1 E(k)) 



1 + A 2 V\ sin 2 kj + (AE(fe) 



1) 



2(AS(fe) + 1) 



E(fc) = m + (1 - cos kj) 



(13) 

(14) 
(15) 



An useful quantity is the mechanical compressibility that, at zero-temperature coincides with the density 
of states: 



1 



d 2 s 

(2tt) 3 ||V fe £(fc)|| 



(16) 



The density of states of the system present a typical band structure (see the upper part of Fig. [j], 
dashed line) . The upper limit of the band corresponds to the saturation due to the Fermi statistics (one 
particle per lattice-site). Since the function E(k) is periodic, its gradient has zeroes in the Brillouin zone, 
producing non-analiticies as the cusps in Fig. (the so-called Van-Hove singularities fl^|). 



3 The Electrical Conductivity 

In a classical paper, Kohn [Q developed an elegant characterization of a conductor, at zero tempera- 
ture. His method allows the measurement of the following limit for the imaginary part of the electrical 
conductivity, a", 

Z = lim Lua"(uj) . (17) 

If this limit turns out to be non zero, the system is a conductor. The construction is as follows. The 
system of interest is constrained to verify periodic boundary conditions in the (say) first spatial direction, 
and immersed in an Aharonov-Bohm like electromagnetic field A = (0, a, 0,0). With this choice of 
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Figure 1: Density of states (upper part) and residue of the conductivity (lower part) for fi > 0. The dashed lines 
correspond to analytical calculations in Fourier space (see Eqs. (|l3|), ([Tq), and (p4|)). The solid lines are obtained 
numerically in finite lattices (see section 0). In all cases r — 1. 



boundary conditions the product L\a a a is gauge invariant since it represents the magnetic flux traversing 
the system. It can be shown that 



Z = 



V s da 2 



(18) 



where Eq is the ground-state energy and V s is the spatial volume. It is crucial that the infinite limit 
volume is taken after the a derivative is performed, since the effect of the Aharonov-Bohm field can be 
thought of as a change in the boundary conditions (see below). In the infinite volume limit, the energy 
no longer depends on a. 

In our case, as the free energy and the ground-state energy coincide in the zero temperature limit, we 
can study the residue in the following way 



lim lim — tr 

Vs^ooT^o dor 



(19) 



where / is obtained from the intensive free energy / after subtracting the vacuum contribution: /(/x) = 

/M-/(o). 

Let us sketch the calculation. The free energy should be calculated in a finite volume and at finite 



4 



temperature. We introduce our system in the Aharonov-Bohm electromagnetic field: 

U wfi = U x<2 = U Xi3 = 1 , U xA = e ia . (20) 
This field can be transformed into a boundary effect by performing the following gauge transformation: 

U x .u - U%„ = e iff(a) U x , v e-^ x+ ^ , g(x) = axt , (21) 

so that U G = 1 excepting 

U$ tl =L 1 -i),i=j aLl - ( 22 ) 
By direct inspection of the fermion matrix in Eq. (^]), one can easily recognize that a system verifying 
periodic boundary conditions in the 1 direction in the field U G is equivalent the same system with no 
field at all, but verifying 

^(x ,x 1 +L 1 ,x 2 ,x 3 ) =e iaLl ¥(x ,xu x 2 ,x 3 ) , (23) 

This amounts to substituting fci by k\ + a in the momentum-quantification in a finite lattice. For a 
system of free fermions the free energy can be now straightforwardly calculated. Once the a derivative 
is performed, the zero temperature limit can be taken by transforming the ko sum into an integral. We 
get, in the simplest case r t = r s — 1, /x > 0, 

Notice that for the empty system, fi < X~ 1 E m i n , the integral vanishes, as well as for the full band 
/i > A _1 i? max , since E(k) is a periodic function of k\. The three dimensional integrals (|24|) can be 
performed using a Monte Carlo method. The results are shown in Fig. |l| (dashed line in the lower part). 



4 Numerical Calculations 

In this section, we are going to reproduce the results of the sections |^ and || by directly considering 
the integration of the partition function. This method has the advantage of being generalizable to 
inhomogeneous external fields, and also when interacting dynamical fields are present. Examples of how 
to introduce an external field on an interacting lattice-gauge system can be found in Refs. Jl3| ]. To 
compute the partition function it is necessary to work in finite lattices, consequently, an infinite volume 
limit should be taken. 

We have carried out measures in symmetric lattices of sizes L = 4, 6, 8, 10, 12, 14 and 16, with m = 1/2 
and A = 1. For the hopping term, we have taken r s = r t = 1. As the integral in the fermionic fields is 
Gaussian, the computation of the electric current just requires the inversion of a AV matrix, V being the 
space-time volume. The fermion matrix (|^) being sparse, we have used a conjugate-gradient algorithm 
for the numerical inversion. 

We first consider the density of states in a vanishing external field. In order to measure dp/dfj, we 
invert the matrix at /i ± e for e small enough. In interacting systems the derivative can be calculated in 
terms of connected correlation- functions [|l Of] . 

The numerical results are plotted in Fig. |l|, upper part, together with the infinite volume values 
obtained analytically. Although the finite size effects are non negligible even in the larger lattices for 
most values of [i, there is a clear trend to the asymptotic values. 

Unfortunately, for an interacting system it is not immediate how to implement Kohn's method for 
calculating the residue of the conductivity. In fact, the free energy is rather hard to calculate with a Monte 
Carlo simulation and what one directly obtains are mean-values. We are now going to present a different 
way of computing the residue, by directly measuring the system response to an external electrical field. 
Notice that the presence of an electric field requires a non-homogeneous vector potential and consequently 
the inversion of the fermion matrix can no longer be performed in closed analytical form. This new recipe 
can be straightforwardly generalized to interacting systems, but its equivalence with the Kohn's method 
is just an ansatz. Nevertheless the agreement is excellent, as we will show. 
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By analogy with continuum electrodynamics, we want to study the electric current induced in the 
system by an external weak uniform electric field in the 1 direction. The conductivity (in the frequency 
domain), will be the proportionality constant between the electrical current and the external field. 

There are some subtleties that need to be considered when putting an external electric field on the 
lattice. We take the gauge-field configuration (t = xq) 

U x , = e i£tx \ U x>i = l, (25) 

£t = T nt, ,-— + 1,..., — -1,— | . (26) 

Notice that the quantization of the electric-field is due to the spatial boundary conditions. To preserve 
the translational symmetry, the displaced gauge field 

U x , = e is ^-$ , £ integer, (27) 



should be a gauge-transform of the one in Eq. (25). Since the needed gauge transformation is analogous 
to Eq. ( pl| ) , it is easy to check that the condition that allows this transformation is the trivialness of the 
Polyakov loop: 

L -l Lq-1 

II C/ (*,-),o = 1 OT Yl £ t = l™, (28) 

t=0 t=0 

with n integer. This condition also allows to transform the gauge field to the Coulomb gauge Aq = 0. If 
condition ( p8|) is violated, the translational invariance is lost and the electric current is no longer spatially 



homogeneous even on a homogeneous electric-field. However, with the correct field choice (28), we get a 
homogeneous electrical current aligned with the external electrical field, and imaginary as anticipated in 
Eq. (|TJ). In order to directly compare with the results obtained with Eq. (p3), let us define 

3(t) = i(7x,i> (29) 

If we want to stay within linear- response theory, we have to postulate a linear relation between the Fourier 
transform of the electrical current j (t) and the external electrical field £t : 

j(cj) = <r(w)i(u) . (30) 

Notice that both j(t) and £t being real, er(— u>) = a*(u>). However, the results can be more cleanly cast 
in terms of a modified Fourier transform for the electrical field: 

S(u) = -L= V Ste-^+W . (31) 

The rationale for this is that the electrical field £t on the lattice lives mid-way between sites at times t 
and t + 1. The modified conductivity a(uj) = ](lu)/£(uj) is related with the previous one by 

cr(w) = a{u)e lw/2 . (32) 

The nice feature of cr(w) is that it turns out to be purely imaginary. 

In Fig. ^ we plot the imaginary part of a(ui) as obtained from a field with hq = I and n\ = — 1, (from 
now on we shall only indicate the non- vanishing rij's), in a system of Wilson fermions with m = 1/2, 
?' = I and jj, = I, that is within the band energy-range and therefore with a non- vanishing Fermi surface. 
We see that for large frequencies the thermodynamic limit is reached in rather small lattices. However, 
at the minimal reachable frequency (2tt/Lq) the conductivity is rapidly growing suggesting a singularity 
at zero frequency. In fact, for a (classical) system of free particles of density n we expect that cr(u>) will 
behave as 

2 

^.free, classical ^ j ^ (33) 

muj 

Notice that if (j(oj) has a pole at w = with a purely imaginary residue the same will hold true for d(u>), 
and both residues will be equal. Although the Euclidean conductivity a(u) do not match the real-time 
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Figure 2: The imaginary part of the conductivity a(uj) of a system of free Wilson fermions at r = 1 , m = 1/2 
and yL = l. 



one (being imaginary, it cannot fulfill the Kramers-Kronig relations), one can formally expect the residues 
to coincide in the passage from u) to iw. This suggest to define the following quantity which will be the 
basic object of our study: 

Z E = luj min a(Lj m[n ) , Lu m[n = — . (34) 

In the Li, Lq — > oo limit, Z E tend to the residue of the pole. In order to measure this, we have considered 
the smallest of possible external disturbances: {no = l,n io / 2 = — I}. 

Our result is shown in Fig. [l], lower part. We see that the Euclidean residue follows quite closely 
Kohn's result, which in fact can be considered as the infinite volume limit for our calculation. Moreover, 
the physical picture is rather transparent: when the band is full, the system gets almost inert, while 
when the band is empty, it can be excited by the external field creating a hole in the Dirac sea. Since the 
smallest possible excitation has frequency 2tt/Lq, to be compared with a gap 2m, it is reasonable that 
at /i = 0, the larger is the space-time lattice, the smaller is the system response. In fact, notice that in 
Fig. [l] when /i is below the lower band limit, the curves get horizontal: in this range of // the system can 
be only excited by crossing the gap between the Dirac sea and the conduction band. And the gap is, of 
course, /x independent in a non-interacting system. 

We remark that our results have been found within the linear response approximation. We can control 
this approximation in several ways. One is to study the Fourier transform of the current, for frequencies 
at which the Fourier transform of the external-field vanishes. In Fig. [| we show the zero-mode of the 
electrical current for the electrical-field {tiq = l,n\ = —1}. We see that this non-linear effect tends to zero 
with growing lattice-size, which is quite reasonable since the minimum possible electric field is 2-k /L\. The 
non-linear corrections are oscillating, but modulated by a rapidly decaying function. Roughly speaking, 
for the largest lattice the non linear effects are of the same order as the distance to the thermodynamic 
limit. A further check can be done by comparing the residue obtained from the data in Fig. |^ with the 
one in Fig. [|: in the L = 14 lattice, the differences are at the 0.3% level, while in the L — 6 lattice the 
differences are at the 1.6% level. Therefore, we believe that non- linear effects are under control for the 
not too small fields that we can deal with. 
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Figure 3: The Fourier transform at zero frequency of the electrical current as a function of the inverse lattice 
size. 

5 Conclusions 

We have presented a simple way of studying the electrical conductivity of a system of Wilson fermions 
at finite density and zero-temperature in a path-integral formalism. In particular, we have computed the 
residue of the zero frequency pole of the conductivity, by numerically considering the linear response to 
an external electric field, varying in Euclidean time. The results have been contrasted with an analytical 
computation based on a method proposed by Kohn, and an excellent agreement has been found. As a 
further cross-check, we have computed the density of states both analytically and numerically in a finite 
lattice, obtaining a nice thermodynamic limit convergence. It should be emphasized that in contrast with 
the analytical calculation which can only be done for a non interacting system (or, at most, for simple 
external fields), the numerical calculations are easily generalizable to more complex models, as fermions 
self-coupled with quartic interactions or via a dynamic bosonic field. 

An open, very interesting question is the possibility of extracting the full real-time electrical conduc- 
tivity function from its Euclidean counterpart. We have shown that the residue of the zero- frequency 
pole can indeed be obtained. It would be also very interesting to extend this approach to systems at 
finite-temperature. 
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